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construct a semiclassical approximation classify themselves into either topo- 



• logical (instantons) and non-topological (bounces) solutions. The quantum 

■ amplitudes consist on an exponential associated with the classical contri- 
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qq ■ bution multiplied by the fluctuation factor which is given by a functional 
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determinant. The eigenfunctions as well as the energy eigenvalues of the 

■ quadratic operators at issue can be written in closed form due to the shape- 

Q 

invariance property. Accordingly we resort to the zeta-function method to 
compute the functional determinants in a systematic way. The effect of the 

• multi-instantons configurations is also carefully considered. To illustrate the 



instanton calculus in a relevant model we go to the double-well potential. The 



second popular case is the periodic-potential where the initial levels split into 
bands. The quantum decay rate of the metastable states in a cubic model is 
evaluated by means of the bounce-like solution. 
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I. INTRODUCTION. 



The tunneling through classically forbidden regions represents one of the most striking 
phenomenon in quantum theory and therefore plays a central role in many areas of modern 
physics. On the other hand, together with the operator formalism of quantum mechanics we 
have an equivalent description by means of path-integrals. In such a case the Schrodinger's 
equation is substituted by a global approach where the quantum mechanical time evolution 
is analysed in terms of a functional integration. Qualitatively speaking, the path-integral 
representation corresponds to a sum over all histories allowed to the physical system we 
are dealing with. To be precise, we need to take into account an imaginary exponential of 
the classical action associated with every path which fulfill the appropiate initial and end 
points conditions. Of course, the quantum amplitude so built is difficult to handle due to 
the oscillating character of the exponential at issue. To avoid problems of this sort we carry 
out the change t — > — ir (known in the literature as the Wick rotation). In doing so, we can 
take advantage of the euclidean version of the path-integral which represents by itself a new 
tool for describing relevant aspects of the quantum theory. 

Almost from the very beginning of the subject a semiclassical treatment of the tunneling 
phenomena (ranging from periodic-potentials in quantum mechanics to Yang-Mills models 
in field theories) has been performed by means of the so-called instantons. Going to more 
physical terms, the instantons represent localised finite-action solutions of the euclidean 
equation of motion. To be specific, the euclidean equation is the same as the usual one for our 
particle in real time except that the potential is turned upside down. Although more massive 
than the perturbative excitations, the instantons themselves become stable since an infinite 
barrier separates them from the ordinary sector of the model. The stability is reinforced by 
the existence of a topological conserved charge which does not arise by Noether's theorem 
in terms of a well-behaved symmetry, but characterizes the global behaviour of the system 
when the imaginary time is large enough. Accordingly, it comes by no surprise that these 
classical solutions have been considered in the literature under the name of topological 
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configurations. Once the appropiate classical solution is well-known, we make an expansion 
around the topological background to evaluate the quadratic fluctuations which arise in 
terms of the functional determinant of a second order differential operator. The integration 
is solved within the gaussian scheme except for the zero-modes which appear due to the 
invariances of the model. To deal with these excitations we introduce collective coordinates 
so that ultimately the gaussian integration is carried out only along the directions orthogonal 
to the zero-modes. 

As a functional determinant includes an infinite product of eigenvalues, the result should 
be in principle a highly divergent expression. Fortunately we can regularize the fluctuation 
factors by means of the ratio of determinants. The eigenfunctions as well as the energy 
eigenvalues of the quadratic operators which appear in this article are obtained in closed form 
by virtue of the shape-invariance property just derived in one-dimensional supersymmetric 
quantum mechanics. Once we reach this point, the evaluation of the quotient of determinants 
is borne out by the zeta-function method. 

On the other hand the use of metastable states, defined as long-lived but eventually 
decaying states, is a fact of crucial relevance in the most diverse branches of physics. The 
classically stable configurations which allows the system to decay through quantum tunneling 
into the true vacuum receives in the literature the name of false vacuum. 

For instance, particular interest has been devoted in recent years to the study of cos- 
mological models where the metastability properties are of great importance to explain the 
different steps of the evolution of the universe as a whole. As regards the very early stages of 
our world, it is customarily assumed the existence of an inflationary era in which the energy 
density is highly dominated by that of the Higgs field trapped precisely in a false vacuum. 
As expected, the end of the inflation period occurs once the metastable state decays to the 
true ground-state of the model. 

Restricting ourselves to quantum mechanics, let us consider a particle moving under the 
action of a potential which exhibits a metastable well around the origin. In such a case, the 
particle will stay there forever as far as classical mechanics concerns. Because of tunneling 
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however, the potential allows the escape of the particle from the single well to infinity. In 
other words, the state situated around the origin is quantum mechanically unstable and 
the probability current outside the barrier is non-vanishing. Taking into the account this 
fact, the energy eigenvalues become complex so that the imaginary part itself represents a 
measure of the decay width. 

In order to compute the lifetime of the particle trapped within the single trough of the 
potential different methods have been proposed. Again the semiclassical expansion based 
on the euclidean version of the path-integral represents the most suitable approach in this 
context. When constructing a semiclassical theory for the decay of metastable states, the 
euclidean non-topological configuration which leads the process is often referred to as a 
bounce solution. To be precise, the spectrum of small fluctuations around the bounce itself 
contains a negative eigenvalue which becomes ultimately responsible for the metastability 
of the system as a whole. In summary, both topological (instantons) and non-topological 
(bounces) euclidean configurations represent probably the most adequate tools to study 
tunneling phenomena. 

It remains to identify the general class of potentials for which the functional determinant 
can be computed according to the method based on the zeta-function. Our description 
takes over the shape-invariance properties of the one- dimensional Schrodinger equation. 
Among other things, the shape-invariance represents a sufficient condition for the solvability 
of the Schrodinger equation at issue 0. In such a case the energy eigenvalues as well as 
the eigenf unctions can be obtained in closed form so that the explicit evaluation of the 
associated zeta-function is feasible. The curious reader can find in [0J a complete list of 
solvable potentials giving rise to shape-invariance properties. To be more specific, the cases 
considered in this article belong to a general class of Posch- Teller potentials labeled by I 
(£ = 1,2,...). For the sake of shortness only the first members of the series are discussed in 
detail although the potentials associated with higher values of I are relevant from a physical 
point of view ||. 

On the other hand, it may be interesting at this point to compare the results obtained 
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according to the one-loop approximation of the path integral and the ones derived from the 
standard WKB method. Almost from the very beginning of the subject the discussions have 
been based on the double-well potential and the periodic sine-Gordon model. However, the 
exact spectrum of the Lame potential has been found recently ||. Although it is not apparent 
a priori how to handle more complicated situations we can conclude that the aforementioned 
case represents an excellent benchmark to test the standard WKB procedure. In order not to 
clutter the paper we refer to the article of Miiller-Kirsten et al. {§] where the interested reader 
could find a careful discussion. For background we simply point out that only the so-called 
WKB-related method of matched asymptotic expansions yields the exact instanton results. 
In other words, we need to go beyond the well-known linear connection relations to reach the 
correct values. To sum up, it may be worth spelling out how the one-loop approximation 
based on the instantonic approach provides a relevant framework for evaluating physical 
magnitudes in a more reliable context than the simple WKB method. 

The arrangement of the article is as follows. First of all, we review in detail the main 
features concerning the instantonic approach in one-dimensional particle mechanics. To 
illustrate the method we resort to the double-well potential. In addition, we also study 
the harmonic oscillator since represents the reference for all the models discussed along 
our analysis. The quantum effects of the one-instanton as well as the multi-instantons 
configurations are included in this chapter. 

The following section is devoted to the periodic-potential. Among other things, the 
lowest band arising from the splitting of the initial ground-state is explained by means of 
the multi-instantons contribution. We conclude with the consideration of a cubic model to 
stick out the existence of metastable states. Different mathematical results concerning one- 
dimensional supersymmetric quantum mechanics, shape-invariance properties, zeta-function 
regularization procedures and spectral densities can be consulted in the appendices. 
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II. THE INSTANTONIC APPROACH. 



To start from scratch let us describe in detail the instanton calculus for the one- 
dimensional spinless particle as can be found for instance in ||. The interested reader 
finds there a comprehensive description of the whole subject. In order to be specific, we 
assume that the particle moves under the action of a confining potential V(x) which yields 
at quantum level a pure discrete spectrum of energy eigenvalues. Unless otherwise noted we 
choose the origin of the energy so that the minima of the potential satisfy V(x) = 0. If we 
set the mass of the particle equal to unity for notational simplicity, the lagrangian L which 
governs the behaviour of the model is given by 

H(f) 2 - 1 ^ « 

Now we can start to quantize the theory. If the particle is located at the initial time 
ti = —T/2 at the point Xi while one finds it when tf = T/2 at the point Xf, the functional 
version of the non-relativistic quantum mechanics allows us to express the transition ampli- 
tude in terms of a sum over all paths joining the world points with coordinates (—T/2, xi) 
and (T/2, xf). At this point it proves convenient to write the action S starting from the 
lagrangian L, i.e. 

S = / L(x,x) dt (2) 

J-T/2 

so that the contribution of a path has a phase proportional to the action S itself. Therefore 
we reduce our problem to the study of a transition amplitude expressed as 

< x f \ exp(-iHT)\ Xl >= N(T) J [dx] expiS[x(t)] (3) 

It may be worth spelling out that H represents the hamiltonian of the model at issue 
while the symbol [dx] indicates the integration over all functions which fulfill the adequate 
boundary conditions. The factor N(T) will be chosen to normalize the amplitude conve- 
niently when we discuss more mathematically the meaning of (^]). As the hamiltonian H 
gives rise to a pure discrete spectrum of energy eigenvalues, namely 

6 



H\n >= E„\n > 



(4) 



we can write that 



< Xf \ exp(—iHT)\xi >= exp(— iE n T) < Xf\n >< n\xi > 



(5) 



As we are mainly interested in the first eigenfunctions of H, it proves convenient to 
transform the exponentials of (§) into decreasing exponentials. For such a purpose, we make 
the change t — > — it, known in the literature as the Wick rotation, so that 



iS[x(t)] 



T/2 
-T/2 



1 / dx s 



- V(x) 



2 \dr 

In the following we resort systematically to the euclidean action S e , i.e. 



(6) 



S P 



T/2 



1 / dx* 

2 ['dr 



V(x) 



dr 



(7) 



To sum up, the euclidean formulation of the path-integral corresponds to 
< Xf \ exp(—HT)\xi >= N(T) J [dx] exp —S e [x(r)] 



(8) 



To expose the main features of the semiclassical approximation we start by considering a 
particular trajectory x c (t) which satisfies the boundary conditions at issue. Next we perform 
the expansion of a generic x(t) with identical Dirichlet conditions as x c (t) according to 



XT = X, 



(9) 



where Xj{r) stand for a complete set of orthonormal functions, i.e. 

r T/2 



T/2 



Xj(r) x k {r) dr = 5 jk 



(10) 



vanishing at our boundary 



Xj {±T/2) = 



(11) 



We carry things further and make the choice 
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[<fe]=n4r (12) 

so that as a matter of fact each path is completely characterized by the Cj themselves. It 
seems plausible to explain the functional formalism in terms of the integration over the 
Fourier coefficients Cj. Now we need to identify the quadratic differential operator which 
gives rise to the basis Xj{r). 

The quasiclassical approximation (or steepest descent method in more mathematical 
language) takes for granted that x c (t) represents a stationary point of the euclidean action. 
As corresponds to an extremal path x c (t) verifies the equation 

^ = V ' {X) (13) 

where the prime denotes as usual the derivative with respect to the spatial coordinate. 
Notice that ([i~3| ) corresponds to the euclidean equation of motion for the particle once the 
potential has been turned upside down. To take into account the quantum fluctuations we 
perform a functional expansion about x c (t). As expected the crucial feature is the analysis 
of the second variational derivative since the linear term is absent since x c (t) represents an 
extremal path. We find that 



rT/2 

S e [x c (r) + 5x(t)} =S eo + / Sx 

J-T/2 



.^ Sx+ i v » lXc{T)]Sx 



dr (14) 



being S eo the classical action associated with the configuration x c (t). The conventional form 
of the semiclassical approximation takes over a complete set of eigenfunctions (eigenvalues) 
of the so-called stability equation, i.e. 

d 2 



dr 2 



v J (r) + V"[x c (r)]v 3 (r) = e J v J (r) (15) 



Now the expansion of (|14]) becomes diagonal so that 



Se = S eo + \ C fa ( 16 ) 



2 

3 

and therefore the euclidean transition amplitude reduces itself to 
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<x f \ exp{-HT)\ Xi >= N(T) exp(-S eo ) Re] 



1/2 



(17) 



Almost from the very beginning of the subject the product of eigenvalues has been written 



as 



n 



.-1/2 



Det 



d 2 



dr- 



+ V"[x c {t)\ 



-1/2 



(18) 



according to a notation which obviously has its origin in the finite-dimensional case. In 
principle we assume that (|Tj|) represents a formal expression where all the eigenvalues €j 
are positive. The difficulties associated with the existence of zero-modes can be avoided by 
introducing collective coordinates, not to mention the negative eigenvalue which yields the 
metastability phenomenon (more on this later). 

Now we remind the properties of the euclidean transition amplitude for the well-grounded 
problem of the harmonic oscillator since this model is the reference to deal with the quotient 
of functional determinants when we go to more relevant systems. If we are dealing with 



V(x) 



v 



X 



(19) 



the amplitude between Xi = and Xf = is given by 

<x f = 0\ exp(-H ho T)\ Xi = >= N(T) ^Det 



dr 2 



-1/2 



The interested reader can find in || the final form of the amplitude, namely 
<x f = 0\exp(-H ho T)\xi = >= Q ' (2 sinh z/TT 1/2 



(20) 



(21) 



Next we must carry things further in order to discuss in detail models which give rise to 
relevant instantonic configurations far beyond the trivial cases where x c {r) reduces itself to a 
constant. (A complete description of the subject can be found in the book of Rajaraman 0). 
To illustrate these ideas in a simple context we again resort to the double- well potential given 
by 



K(x) = ^(x 2 -1) 2 



(22) 



The model is particularly well adapted for such a purpose since at classical level we find 
two degenerate minima located at x_ = —1 and x + = 1 as expected when considering a 
potential which enjoys the discrete symmetry x — ► —x. 

The question we wish to address now is the explicit description of the tunneling phe- 
nomenon in the euclidean version of the path-integral. To make clear the method we turn 
upside down the potential as required by the Wick rotation. From a physical point of view 
it seems reasonable to evaluate the transition amplitude between the points x_ = — 1 and 
x+ = 1. As regards the topological configuration which represents the basis for a quasi- 
classical approximation, we look for a trajectory with initial Xi — — 1 at ij = — T/2 while 
Xf = 1 when tf = T/2. In other words, we need to know in closed form the instanton just 
interpolating between the two minima of the potential V(x). It is customarily assumed that 
T — > oo mainly because the explicit solution of the problem is much more complicated for 
finite T (more on this point later). 

With the appropiate boundary conditions we must integrate (p~3|) as far as the double- 
well concerns. The x c (r) corresponds to the motion of a particle whose conserved euclidean 
energy E is given by 

2 V dr 



E=U d -^) ~V(x e ) (23) 



According to the Dirichlet conditions required by the form of our transition amplitude 
we conclude that E = 0. As a matter of fact the particle leaves the point X, = — 1 and takes 
an infinite time to climb up to the top of the mountain Xf = 1. Once the journey starts 
the particle cannot return to the initial point or explore unlimited motions towards plus or 
minus infinity. The only dynamical possibility is to reach Xf — 1 asymptotically as tf — > oo. 

The aforementioned mechanical analogy allows us to solve the problem just by integra- 
tion of a first-order differential equation instead of resorting to the conventional euclidean 



equation of motion written in (|13[). Therefore we have 

dx 



(h ±pV(x c ) (24) 
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where we recognize the mechanical version of the well-grounded Bogomol'nyi condition first 
introduced in the mid-seventies for the analysis of field theories |8| . Now we explicitly solve 
(p^) by a simple quadrature so that 

lo(t - T c ) 



x c {t) = tanh ■ 



(25) 



where the parameter r c indicates the point at which the instanton makes the jump. As 
expected the antiinstanton is obtained by means of the transformation r — > — r and allows 
the connection between = 1 and Xf = —1. Now we must evaluate the action associated 
with x c {t) according to (0) to get 



2lo 

T 



.V,„ = — (26) 
It should be emphasized that in principle we need classical configurations for which 



x = ±1 at large but finite values r = ±T/2. However the instantons that appear in the 
literature correspond to solutions where the particle takes itself an infinite time to complete 
the trajectory. Fortunately the difference is so small that can be ignored mainly because 
ultimately we are precisely interested in the asymptotic limit T — > oo. As corresponds to a 
semiclassical approximation we now evaluate the quadratic corrections over the background 
provided by the classical configuration. In other words, we take into account the functional 
determinant at issue while the cubic and higher terms in the quantum fluctuations are 
neglected. First of all we need an operator of reference for the explicit evaluation of the 
determinant. The most natural choice is of course the harmonic oscillator just located at 
either of the two minima of V(x), namely 



v"(x = ±r 



OJ 



(27) 



so that our description takes over 



<x f = l\ exp(-HT)\xi = -1 >= N(T) ^Det 



dr A 



-1/2 



Det [-(d 2 /dT 2 ) + V"[x c (r)}} 



-1/2 



Det [-(d 2 /dr 2 



exp(-S' e 



(28) 
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where we have multiplied and divided by the determinant of the aforementioned harmonic 
oscillator. Incorporating the explicit form of x c (r) the stability equation is given by 

d 2 



dr 2 



Xj(t) + 



UJ 



2 3^ 2 



xj(t) = e jXj (T) (29) 



2cosh 2 (cjr/2)_ 

and corresponds to a Schrodinger's equation with exactly solvable Posch- Teller potential. A 
complete description about the eigenfunctions and the eigenvalues of (p9|) is possible in terms 
of the shape- invariance properties which appear in supersymmetric quantum mechanics [|TJ. 
This modern framework is particularly well adapted to incorporate the zeta-function method 
||. (The technical details of the method can be found in appendix A). 

Among other things one finds a zero-mode x q (t) which could jeopardize the physical 
meaning of the determinant. However this eigenvalue e G = comes by no surprise since it 
reflects the translational invariance of the system. To be more specific, there is one direction 
in the functional space of the second variations which is incapable of changing the action. 
One could just as well discover the existence of a zero- mode starting from (0). If we perform 
an additional derivative with respect to r, then 

= VCJF (30) 



is just the solution of (jig ) with e Q = 0. Notice the normalization of x (r) which is due 
precisely to the zero euclidean energy condition for x c (r). Going back to the functional 
arguments we can see how the integration over the Fourier coefficient c G becomes tantamount 
to the integration which takes care of the center of the instanton r c . To fix the jacobian of 
the transformation we consider a first change like 

Ax(r) = x {t) Ac (31) 

According to the general expression written in (Q) we observe that under a shift Ar c the 
effect corresponds to 

Ax(r) = -JsZ x {t)At c (32) 
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Next the identification between and fl32|) yields 

dc Q = \fsZdT c (33) 

where the minus sign disappears since what matters is the modulus of the jacobian at issue. 
In doing so we have that 

\ Det [-(d 2 /dT 2 ) + V"[x c (r)}} }- 1/2 _ 
\ Det[-{d 2 /dr 2 )+uj 2 } J 

i Det'[-(d 2 /dT 2 ) + V"[x c (r)]] \- 1/2 fsZ ( ^ 
\ Det[-(d 2 /dr 2 )+^] J V 2tt c 1 > 

where Det' stands for the so-called reduced determinant once the zero-mode has been ex- 
plicitly removed. 

Different procedures have been used in the literature to evaluate the typical ratio of 
determinants which appears in fl3"l]). The reasonings become much more accesible once 
the system is enclosed in a box of lenght T, thus avoiding the subleties associated with 
the continuous spectrum. If one imposes periodic boundary conditions the most relevant 
physical information derives from the analysis of the phase shifts [f| . At last the limit T — > oo 
is necessary to achieve a meaningful result. However we prefer from the very beginning to 
work in open space and compute the aforementioned ratio of determinants by resorting to 
the zeta-function regularization procedure [0. Going back to fl29|) and (0), we perform 
the change of variable 

; = ^ (35) 

to recognize the presence of the ratio of determinants 

n Det ' ° 2 

Q2 = ~D^T 2 (36) 

together with a factor (3 given by 

P = £ (37) 
13 



As regards the non-zero spectrum of the operator O2 we find a discrete level with Ei = 3 
while the energy of the scattering states corresponds to 



E k = k 2 + 4 



(38) 



for eigenfunctions derived from the plane waves according to 



<p2,k{ Z ) 



4(z) A\(i 



exp(ikz) 



2n 



If we bear in mind the form of the operators A\(z) and A{(z), namely 

A\ — — ; — h tanh z 

dz 



(39) 



(40) 



Al 



dz 



2 tanh z 



we can write that 



-k + 2 



3ik tanh z 



exp(ikz) 



(41) 



(42) 



cosh^z J V2n y/k 2 + 4 yjk 2 + 1 

As a straightforward consequence of fl42|) the regularized spectral density p r (k) of our 
problem reads 

3(k 2 + 2) 



p r (k) 



(43) 



7T(fc 2 +4)(fc 2 + 1) 

In this scheme the suitable zeta-function ( r {s) to compute the ratio of determinants of 
]) can be expressed as 



Cr(s) = Co 2 ( s ) - (p 2 (s) 



(44) 



Armed with this information we write 



r( , 1 [°° [ . 3 r°° (fc 2 + 2)exp[-(A; 2 + 4)/i] 

US) = fls) Jo " ^ 6Xp( " 3/i) - * loo (fc 2 + 4)(fc 2 + l) 



dk 



(45) 



which, in turn, gives rise to JTT 

Cr(s) 



1 3 /■<» (k 2 + 2) 



3 s ttJ-oo (Jfe 2 + l)(Jfe 2 + 4)»+ 1 
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(46) 



Next it suffices to break the integral of ( fl6| ) into more simple components to obtain ( r {s) 
in terms of Gamma and Hypergeometric Functions, namely |TT[ 

1 3 1 r ( s + l) 3 1 r ( s + |) / 3 3\ 
Cr(s) = 3^ " 7?2^+i r(s + l) " 0F2^+3 T( S + 2) F I 1 ' S + 2' * + 2 ' 4 J (4?) 



If we incorporate that [12 



K44H < 48 > 

the evaluation of £ r (s) at s = yields 

Cr(0) = -1 (49) 

On the other hand the interested reader can find in the appendix B the fundamental 
steps which allow ultimately the computation of Next if we simply collect the partial 

results obtained along the previous paragraphs the ratio of determinants R which appears 
in ©, i.e. 

Det'[-(dyd^) + V"[x c (r)}} 

Det[~(d 2 /dr 2 )+uj 2 } 1 1 

reduces itself to 

R= ih (51) 

Combining now these informations we obtain the transition amplitude anticipated in 
(f28|), namely 

< Xf = 1| exp(— HT)\xi = —1 >= 

(2sinhu/rr 1/2 ^J~S^ J^- exp(-S eo ) u dr c (52) 

Apart from the first factor, which represents the contribution of the harmonic oscillator 
of reference, we get a transition amplitude just depending on the point r c at which the 
instanton makes precisely the jump. According to the values of the interval T the result 
seems plausible whenever 
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'S eo y ^ exp(-5' eo ) u; T < 1 (53) 

a nonsense condition if T is large enough. However in this regime we can accommodate 
configurations constructed of instantons and antiinstantons which mimic the behaviour of 
a trajectory strictly derived from the euclidean equation of motion. In doing so we get 
an additional bonus since the integration over the centers of the string of instantons and 
antiinstantons is carried out in a systematic way. As a matter of fact the well-known formula 
for the level splitting of the double-well potential appears when considering the effect of the 
multi-instantons. 

As all the above calculations were carried out over a single instanton, it remains to 
identify the contributions which take into account the effect of a string of widely separated 
instantons and antiinstantons along the r axis. It is customarily assumed that these combi- 
nations of topological solutions represent no strong deviations of the trajectories just derived 
from the euclidean equation of motion without any kind of approximation. We shall com- 
pute the functional integral by summing over all such configurations, with j instantons and 
antiinstantons centered at points t±, ...,Tj whenever 



T T 
-- < n < ... <Tj<- (54) 



Being narrow enough the regions where the instantons (antiinstantons) make the jump, 
the action of the proposed path is almost extremal. We can carry things further and assume 
that the action of the string of instantons and antiinstantons is given by the sum of the j 
individual actions. This scheme is well-known in the literature where it appears with the 



name of dilute gas approximation [13 



Now we can evaluate transition amplitudes with closed paths with Xi = —1 = Xf for 
instance, so that the action at issue St will be an even multiple of the single instanton 
action, i.e. 

S t = 2] S eo (55) 



16 



As expected the amplitude between Xi = —1 and Xi = 1 incorporates a contribution 
given by 



S t = (2j + 1) S e 



(56) 



In addition the translational degrees of freedom of the separated j instantons and anti- 
instantons yield an integral of the form 



T/2 fTj p (uTy 

-T/2 7-T/2 ./-T/2 j! 



(57) 



As regards the quadratic fluctuations around the j topological solutions we have now 



that the single ratio of determinants transforms into [[13 



7T 



(2 sinhu/T)~ 



1/2 f ^ef [-(d 2 /dT 2 ) + y"[a c (T)]] 
\ Det[-(d 2 /dr 2 ) +u 2 } 



-1/2 



7T 



1/2 



exp(— u/T/2) 



Det' [-(d 2 /dT 2 ) + V"[x c (r)}} 



-1/2' 



(58) 



£>ei[-(d 2 /dT 2 )+o; 2 ] J 

according to the limit of the factor associated with the harmonic oscillator when T is large. 
Next we can write the complete transition amplitudes for the double-well potential so that 



< Xf — 1| exp(— HT)\xi = —1 >= 



exp(-u/T/2) ^ 



7T 



t (2J + 1)! 



where cf stands for the so-called instanton density, i.e. 



d = \l - y S eo exp(—S e 



In summary 



< X/ = 1| exp(— HT)\xi = —1 >= f — J exp(— luT/2) smh.(uTd) 



Similarly 



< Xf = 1| exp(— i?T) |xj = 1 >- 



exp(-a/T/2) ^ 



7T 



and therefore 



(59) 



(60) 



(61) 



(62) 
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< Xf = 1\ exp(—HT)\xi = 1 > 



( 




exp(— uT/2) cosh(uTd) 



(63) 



Resorting to the limit T — > oo in (|6T|), we obtain the energy eigenvalues E Q and Ei of 



the first two levels of the double-well potential 



E, 



2 



2lu J — exp(— 2u;/3) 



(64) 




(65) 



As expected the quantum mechanical tunneling transfers the wave function from one well 
to the other, thus lifting the degeneracy of the classical vacua. To close, notice the way in 
which the energy eigenvalues depend on the barrier-penetration factor, i.e. the exponential 
of minus the classical action of the instanton at issue. Once we have understood the main 
properties of the instanton calculus for the one-dimensional particle we proceed to extend the 
method to well different models like a periodic-potential based on the sine-Gordon theory or 
the cubic system which represents by itself an excellent benchmark to discuss the existence 
of metastable states. 



In this section we consider the particle under the effects of a periodic-potential V(x). At 
first one can ignore the barrier-penetration so that the energy eigenkets are infinitely degen- 
erate. Going to more physical terms, we construct a set of harmonic oscillators centered at 
the bottom of each well. However the quantum-mechanical tunneling changes dramatically 
this naif picture. As a matter of fact the single energy eigenvalue transforms into a continu- 
ous band while the eigenstates are given in terms of Bloch waves. Our objective is to describe 
how the instanton method serves to explain these results which are obtained of course by 
means of conventional procedures in any solid-state physics course. To be more specific, the 
periodic-potential is the most used model for the study of electrons in a one-dimensional 
lattice. 



III. THE PERIODIC-POTENTIAL. 
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Once we have carried out the discussion of the one-instanton amplitude, the dilute-gas 
approximation takes over a set of instantons and antiinstantons freely distributed along 
the r-axis. As expected each topological configuration (instanton or antiinstanton) starts 
just at the point where its predecessor ends. In addition, the number of instantons minus 
the number of antiinstantons represents the observed change in the coordinate x between 
the initial and final points of the transition amplitude at issue. For reference we choose a 
concrete periodic-potential V(x) given by 



1 — cos 



(66) 



V(X) = LU 4 

so that the minima of V(x) lie at 

x = 2nju), j e Z (67) 

Notice that as usual the minima satisfy V(x) = 0. This periodic-potential can be un- 
derstood as the quantum-mechanical version of the well-grounded sine-Gordon model which 
consists of a real scalar field <fi(x, t) in d — 1 + 1 dimensions governed by a potential function 
similar to flBDP . As a matter of fact the instanton we are interested in represents the soliton 
in the corresponding bidimensional field theory. 

From a qualitative perspective the starting point should be the so-called tight-binding 
approximation where the tunneling of the low-energy eigenkets from one well to the next one 
is irrelevant. In doing so we have an infinitely degenerate ground-state with energy given by 

Eo = ^M (68) 

as corresponds to a set of harmonic oscillators distributed along the one- dimensional lattice. 
When considering the tunneling effects the single level E yields a complete band, i.e. 

E e = E a -a cos (69) 

being a a constant to determine while 9 serves to denote the different states along the band 
itself H . 
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Next we explain the way in which the instanton calculus provides an accurate description 
of this Bloch wave. To start from scratch we take into account the transition amplitude with 
Xi = at ti = —T/2 and Xf — 2-kuj when tf — T/2. In doing so the functional integral 
is dominated by the instanton which interpolates between the two adjoining minima of the 
potential V(x). The explicit form of the topological configuration derives once more from 
the Bogomol'nyi condition. Solving (|2~4[) by means of a quadrature we have 



x c (t) = 4co>arctan [exptu(r — r c )] (70) 
so that the action associated with x c (t) reads 

S eo = 8iU 3 (71) 

Next we must write the stability equation over the instanton, i.e. 
d 2 



2 2w 2 



Xj{t) = ejXjij) (72) 



COsh 2 LUT 

to obtain again a Posch- Teller potential. In any case we conclude that 

< x f = 2nuj\exp(-HT)\x i = >= I- J (2 sinhwT) _1/2 



l Det'[-(d 2 /dr 2 ) + V"[x c (T)]] \- 1/2 fsZ m , 
\ Det[-{d?/dT 2 )+u 2 } J V 2tt c { ' 

since V"{2Ttju) = uj 2 while r c represents as usual the collective coordinate of the problem. 
Now we perform the change of variable 

z = cor (74) 
so that the relevant quotient of determinants corresponds to 

n Det ' 0l frt\ 

Ql = doK (75) 

where the global factor is now (3 = uj 2 . The discrete spectrum reduces itself to the zero- mode 
while the scattering states have energies of the form 
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E k = k 2 + 1 



(76) 



for eigenkets like 



4>i,k( z ) 



In other words 



h,k( z ) — (~ik + tanhz) 



exp(ikz) 



exp(ifcz) 



(77) 



As regards the expression of the regularized spectral density p r {k) we find that 

p r (k) - 



(78) 



(79) 



7r(k 2 + 1) 

The structure that is at work corresponds therefore to the zeta-function Cr( s ) given by 



Us) = C 0l (s) ~ CpAs) 



(80) 



To be more specific 



Cr(s) = T^T / ^ ^ — / 
1 (SI JO 7T J- 



1 f°° exp[-(£; 2 + l)fj] 



{k 2 + 1) 



dk 



(81) 



In accordance with (81) we have that 



Us) _M 



r(s + i) 

Taking advantage of the well-known properties of the Gamma Function we obtain 



(82) 



Cr(0)=-1 



(83) 



together with (see appendix B) 



C(0) = 2 In 2 



(84) 



This approach provides the conventional ratio of determinants, i.e. 

1 



R 



4u 2 
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(85) 



which in turn gives rise to the transition amplitude we are looking for. To sum up 

< xj = 2nui\ exp(— HT)\xi = >= 



V 2 , / — F) 

1 . , „x -1/2 / ~ / <^ 



2sinhu;T) ' J S eo \ — exp(—S eo ) u dr c (86) 

Now, when going to the dilute-gas approximation we must consider a set of instantons 
and antiinstantons so that each topological configuration starts where its predecessor ends. 
As regards the transition amplitude fl86|) we require that the number of instantons (n) minus 
the number of antiinstantons (n) fulfill the condition n — n = 1 as corresponds to the jump 
between Xi = and Xf — 2-kuj. Because of indistinguishability the dilute-gas approximation 
includes a combinatorial factor F in order to eliminate the over-counting of configurations 
which correspond to the exchange of center-positions. For the amplitude (0) we find that 

(n + n)\ 



(87) 



Notice the difference with the double-well potential where the instantons strictly alternate 
with the antiinstantons since the problem has only two minima so that the combinatorial 
factor is not necessary at all. 

By integrating over the translational degrees of freedom we obtain now a term I of the 
form 

I = f±l (88) 
(n + n)! 

Armed with this information we can write the transition amplitude in the limit T — > oo, 

i.e. 

<x f = 2ttu\ exp(-HT)\ Xi = >= - exp(-^T/2) ^ 1 1 5 n ^,i (89) 



Try ^ m m 



where 5 n -n,i is the Kronecker symbol while d stands as usual for the so-called instanton 
density, namely 

d = \fs^ exp(-S' eo ) (90) 
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with S eo as given in (|71|) . Now the sum over n and n decouples by resorting to the integral 
representation of the Kronecker symbol, i.e. 

d6 



J n—n,l 



/ — exp[-z0(n - n - 1)1 (91) 
Jo 2n 



so that ultimately we have two independent exponential series. Consequently our main 
result should be 

< Xf = 2ttuj\ exp(—HT)\xi = >= 
/ — exp(z0) - exp(2udT cos 9 - uT/2) (92) 

JO 27T \7Tj 



so that the left-hand side of ( 192 ) is dominated when T — > oo by a continuous band 
parametrised in terms of 9 with < 9 < 2n. Therefore the energy of the states contained 
in this low-lying band should be 

E e = | - 2udcos9 (93) 



As expected the result is consistent with the expression written in ( |69"D so that we have 
described by means of the instanton calculus the behaviour of the lowest band of a periodic- 
potential. 



IV. TUNNELING AND DECAY. 

Next we take into account the decay of metastable states by tunneling processes. In 
principle the states with a finite lifetime are analytically described by means of energy 
eigenvalues which appear in the lower half of a complex plane due to a negative imaginary 
part. According to the temporal evolution of a standard wave function it is the case that 
such imaginary part represents a measure of the decay width F through the relation 

T = -21m E (94) 

As a matter of fact the negative sign of this imaginary contribution for the energy 
eigenvalues serves to prevent the existence of unphysical states with exponentially growing 
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norm. Our purpose along this section should be the computation of the magnitude V for 
the lowest state of a concrete problem. From a physical point of view we assume as usual 
a high barrier-potential so that the lifetime at issue is long. In other words, the state is 
approximate stationary so that a quasiclassical analysis represents the most suitable tool to 
discuss the main characteristics of the system. 

To be more specific, let us consider the point particle in the presence of a cubic potential 

like 

V(x) = yx 2 - l* 3 (95) 

Clearly the potential has a minimum located at x — 0. According to classical mechanics, 
once the particle is trapped in the single trough of V(x) (with negligible kinetic energy) 
the stability is strictly guaranteed. However the situation is dramatically different in the 
quantum regime because of tunneling. As a matter of fact the particle can go towards 
+00 with a non- vanishing probability current there. Consequently non-self adjoint boundary 
conditions must be added to the hamiltonian so that the energy eigenvalues become complex. 
The question we wish to address now is the computation of T for the lowest state of (|95|). 
Qualitatively speaking, the classical configuration which dominates the tunneling process 
is called a bounce. As expected the metastability properties derive from the behaviour 
of the fluctuations built over the classical configuration. When considering non-topological 
bounces the spectrum of the stability equation has one negative eigenvalue which is therefore 
responsible for the decay phenomenon. It may be interesting at this point to remind for 
instance the case of the double- well potential in this context. 

A particle released from rest at the top of the left-hump of —V(x) when ti = —00 arrives 
at the other hump with velocity zero at time tf = 00. In other words, the velocity increases 
from zero to the corresponding maximum and decreases to recover asymptotically the zero 
again. But this velocity represents the zero-mode of the stability equation. Simply put, 
such zero-mode is a nodeless wave function so that the second variational derivative of the 
euclidean action is positive semidefinite. 
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In the semiclassical approximation we shall compute the transition matrix element be- 
tween Xi = and Xf — 0. As usual the development relies on the solution obtained by 
solving the imaginary-time equation of motion which is equivalent to a conventional tra- 
jectory in the reversed potential —V(x). The particle starts at Xi = 0, goes through the 
minimum of —V(x) and then returns to the initial point. Because of this come back motion 
the non-topological solution receives the name of bounce. 

Next we proceed in the same manner as in the previous sections. First of all we resort 
again to the Bogomol'nyi condition to get the form of the bounce-like configuration, i.e. 

*c(t) = 2 Y ~\~~Py\ ( 96 ) 

cosh [lj(t — T c )/2\ 

so that the euclidean action of d96f) should be 



S„ = ^ (97) 

The Schrodinger equation which governs the behaviour of the fluctuations over the 
bounce reads 

d 2 



dr 2 



xAt) + 



2 3lu 2 



X A T ) = W(r) (98) 



cosh 2 ur/2 

In the case of this cubic potential, the velocity of the particle at issue has a zero located 
precisely at the turning point. To sum up, the second variational derivative of the euclidean 
action is expected to possess one lower even eigenket with negative energy eigenvalue. This 
is the point at which the path integral formalism introduces the met ast ability. However the 
explicit derivation of Im E requires a careful procedure. 

First of all we need to perform the integration over the eigenstate x_i(r) with eigenvalue 
e_i < 0. For such a purpose one takes into account the analytic continuation procedure 
through a continuous sequence of paths drawn in the functional space at issue. In order 
not to clutter the paper we omit the mathematical details of the computation to emphasize 
the physical consequences of the result. The interested reader can find in || a complete 
description of the analytic continuation method as tailored to this specific integral. Therefore 
we have that 
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I exp(-e_ 1 c 2 _ 1 /2) 



(99) 



so that the tunneling rate formula will include the presence of the additional factor 1/2 
together with the modulus of the quotient of determinants. According to these arguments 
the transition amplitude over the bounce reads 

'UA 1 / 2 



<x f = 0\ exp(-HT)\xi = 0>- 



7T 



[2 sinhcuT) 



-1/2 



Det' [-{d 2 /dr 2 ) + V"[x c {t)\] 



-1/2 



2tt 



dr r 



(100) 



Det [-(d 2 /dr 2 ) + u 2 ] 

because V"(0) = u 2 . Once we incorporate the change of variable fl35|) the ratio of determi- 
nants now corresponds to 

Det' (0 3 - 5) 



'101 



Det P 2 

with the factor (3 given in fl3~7l) . As regards the relevant discrete spectrum of O3 — 5 we have 
E-i = —5 and E\ = 3 while the energy of the scattering states corresponds again to (p8|). 
As a matter of fact these eigenkets derive from the plane waves, i.e. 



4(z) 4(z) A\(z) 



exp(ikz) 



(102) 



to give 



& 3,fcUJ 



15sinhz ,,, . , / 15 , , 

6(A; 2 - 1) tanhz + ik f — — + k - 11 



cosh z 



Vcosh 2 . 



exp(ifcz) 



Consequently the regularized spectral density p r {k) is of the form 

6(fc 4 + 8fc 2 + 11) 
M } ~ tt(P + 9)(F + 4)(A; 2 + 1) 

As expected the zeta- function ( r (s) reads 

Cr(s) = Co 3 -5(s) - (P 2 (s) 



(103) 



(104) 



(105) 
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so that 



Cr{s) = — - / (jf^dfi 

T(s) Jo 



b f 

exp(— 3/i) + exp(— 5/i) / 

71 J- 

In other words 



6 (A; 4 + 8A; 2 + ll)exp[-(A; 2 + 4)/i] 



(k 2 + 9){k 2 + 4)(P + 1) 



dk 



(106) 



(A; 4 + 8A; 2 + 11) 



- g s + gs - - _/_ oo ^ 2 + g^p + ^2 + A y +1 



dk 



(107) 



Breaking as usual ( |107| ) into more simple integrals we find 
^ r[ > 3* 5* v^2 2 -r( S + l) v^F2 2s +3 T(s + 2) V 2' ' '4/ 



5 1 r(s + §) / 3 5 

+ 7^T^ F ( s+ 2' s + M + 2 '9 



(108) 



If we bear in mind that 



/ ('-,1,2,- 

,2' ' '9 



9 

5 



(109) 



we get again that 



Cr(O) 



(110) 



Taking advantage of the explicit results collected in appendix B the transition amplitude 
at issue reduces itself to 



< x f = 0\exp(-HT)\xi = >= f-J (2smhu/TT 



1/2 



^a/^o \/~~ exp(-S eo ) w dr c 

I s \ 71 



111) 



To close the section we evaluate the transition amplitude by summing over configurations 
which contain a set of separated bounces. In this case we have no restriction to an even or 
odd number of bounces so that the complete exponential series arises. To sum up 
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<x f = Q\ exp(-HT)\xi = >= ( -J exp(-cjT/2) exp(idw/2) 



;ii2) 



for 



'S'eo ex P ( 'S'e 



;ii3) 



In doing so we obtain the final expression of the decay width, i.e. 



'30 



V = u\ — JS eo exp(-S' e 

7T 



;ii4) 



where the prefactor in ( 114 ) is known as the quantum attempt frequency while the second 
term is interpreted as the Boltzman weight for the appearance of the bounce which ultimately 
leads the tunneling phenomenon. 
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V. APPENDIX A. 



In this article we focus the attention on the set of hamiltonians Ot given by 

d 2 £(£ + 1) 



O f 



dz 2 cosh 2 z ~*~ 



(115) 



where £ = 1,2, although only the first members of the series lead to relevant models in 
physics. The operator Oi can be factorized in terms of a superpotential W(z,£) like |j] 



W(z, 



' tanh z 



(116) 



so that Oi = Q\Qi for 



Qi = - — h £ tanh z 

dz 



Qi = — : — h £tanhz 
dz 



117) 



;ii8) 



The shape- invar iance condition appears once we write the partner hamiltonian Oe so 
that in our case the mapping between the old parameter £ and the new one £ reduces to 



- 1 



(119) 



By iteration of the procedure we construct a family of well-behaved superpotentials which 
allow us to solve Oe. The discrete spectrum includes as usual a normalizable zero-energy 
mode (bfJz) of the form 



<t>l,o{ Z ) 



2{2£-l)\ i 



2 e \£-l)\ cosh** 
together with the set of states <f>e >m (z) given by 

1 



(120) 



2(2£-2m- 1)! 
2^{£-m-V)\ J U ^\E m -E : 



Qi(z) . . . Qi_ m+1 (z) 



Lcosh^ m * 



;m) 



with energies 



E m = £ 2 -{£- m f 



122) 
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for m = 1, 



1. The continuous spectrum corresponds to 



H,k(Z 



Q\{z) QU(z) 



Q\{z) 



exp(ikz) 



with energy eigenvalues 



7T 



(123) 



E k = k + 



2 , 1)2 



(124) 



and normalization as 



^lk( z ) <Pe,k>(z) dz = 5{k - k') 



(125) 



Let us now sketch the way in which we can compute the ratio of the determinants 
associated with Oe and Pi which is given by 

d 2 „„ 



dz 2 



(126) 



As a matter of fact Pi itself represents a hamiltonian of comparison for Oe. In such a 
case we need to deal with 

Def Of 



Q 



(127) 



Det P £ 

where Det' Of denotes as usual the reduced determinant once the zero-mode has been 
explicitly removed. To take into account the continuous spectrum we resort to the density 
matrix Hq £ (A;, z, w) written as [|15[ 



E 0e (k,z,w) = <frl k (z) <f> t ,k(w) 



128) 



When going to Pi we get 



T Pe (k,z,w) 



exp —ik(z — w) 
2^ 



(129) 



as expected for a free-particle. The conventional subtraction procedure yields the regularized 
spectral density p r {k) expressed as 



/oo 
[E 0t (k,z,z) - T Pe (k, z, z)} dz 
-oo 



(130) 
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In this scheme the zeta-function for the analysis of ( [[27] ) would be 



Cr(s) = CoA s ) ~ (Pe( s ) 



(131) 



and therefore 



1 r 00 f 
Cr(s) = p^y y o d/i |exp 



E ^ - (< - ™) 2 ]^ 



m=l 



+ / p r (fc)exp -{k 2 + f)fi 



dk 



(132) 
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VI. APPENDIX B. 



Let us come down to the concrete details which allow the evaluation of Cr(0) when 
considering the instanton of the double- well potential. According to the expression of ( r (s) 
(see (Wfl)) we find that 



C(0) = - In 3 + 8 In 2 - \ - ^F'(s + |, 1, s + 2 



\s=0 



(133) 



once we take into account that flT2| 



r'(~) = -0r" ( 7 + 2 In 2) 



(134) 



r'(i) = - 7 



(135) 



r'(-) 



t fl + ln2-l) 



(136) 



r'( 2 ) = -7 + 1 



(137) 



where 7 is the Euler's constant. To obtain the derivative of the hypergeometric function we 
resort to the integral representation, namely |L2| 



v 2' ' '4 ; r 1 r (s + i) Jo v ; v 4 7 



138) 



so that 



F'(s + ^,1, S + 2,^)\ S=0 =h + I 2 + I 3 



(139) 



for integrals of the form 



Jo 



1 3i V 5 , 

1 dt 

4 



(140) 



/ 2 = I'^-j) 5 ln(l-t)cft 



(141) 



32 







If we have that 



4/ V 4 



32 2 
Jo = — m - 
3 3 



16 32, rt 
I 3 = + — In 2 

3 3 3 



then Cr(0) is obviously 



C(0) = 41n2 + ln3 



As regards the bounce solution of the cubic potential we obtain that 
C(0) = 91n2-ln5 + ^ + ^F'( S + | S +l, S + 2^)| s=0 
Taking into account that 

F(s i 3 i.io 5 )- r(s+2) /v d-^ 1 

we find again a situation of the form 



F'(s + -,s + l,s + 2,^)\ s=0 =h + I 2 + I 3 



where in this case we have that 



'-jfO-")" 1 * 



*=jt 1 ( l -!)" ,ht * 



Jo 



5t\~* / 5f\ 
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After a tedious but straightforward computation we can write that 

h = I (153) 

/ 2 = ^ln5--ln2--ln3 (154) 
5 5 5 

54, 54, 18 , , 

h = -— ln2 + — ln3- — (155) 
5 5 5 

In doing so (' r (0) corresponds to 

C(0) = 4 In 2 + In 3 + In 5 (156) 
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